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Abstract 

We present ah initio calculations of the complete gauge-invariant set of two-photon exchange 
graphs for the (ls)^2p3/2 electron configuration in Li-Iike ions. These calculations are an impor- 
tant step towards the precise theoretical determination of the 2p3/2-2s transition energy in the 
framework of QED. 
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I. INTRODUCTION 



At present the lowest-lying states in heavy Li-like ions can be investigated very precisely 
both theoretically and experimentally. One of the most precise experimental results in these 
systems has been obtained by Beiersdorfer and co-workers Q for the 2p3/2-2s transition 
energy in Li-like bismuth, which was determined with an accuracy of 0.04 eV. Accurate 
experimental data are available at present for a number of other elements as well. For 
the latest h,gh-preasio„ measurements we refer to Refs. 333; the outhne of earher 
investigations can be found in Ref. 

The accuracy reached in experimental investigations provides a promising tool for probing 
QED corrections in the strong Coulomb field of the nucleus up to second order in the fine 
structure constant a. For the 2r)i/9-2s transition, this proiect has been carried out in a 
series of our prevous mvest.gat.cr fiflflB. Irr Ref. Q we eomp.eted the eva.uat,o,r of all 
two-electron QED corrections of second order in a and obtained most accurate theoretical 
predictions for the 2pi/2-2s splitting within a wide range of nuclear charge numbers Z. 
Based on a careful estimate of the uncertainty of the theoretical values, we concluded that 
already now the comparison of theory and experiment for Li-like uranium provides a test 
of QED effects of second order in a at the level of accuracy of about 17%. For the 2^3/2- 
2s and 2pi/2-2s transitions in Li-like bismuth, analogous calculations have been performed 
recently by Sapirstein and Cheng However, in order to match with the experimental 
accuracy for the 2p3/2-2s splitting, rigorous evaluations of second-order QED corrections are 
required also for other ions than bismuth. The first step in this direction has been performed 
in our earlier investigation jsl where we have evaluated the vacuum-polarization screening 
correction for several energy levels of Li-like ions, including the (ls)^2p3/2 state. The aim 
of the present work is to calculate the two-photon exchange correction for this state (for 
extensive calculations of these corrections for the lower states in Li-like ions and for non- 
mixed low-lying states in He-like ions we refer the reader to Refs. 3, Is, 13, till, ll^ ll^ Q)- 
After all that, the self-energy screening correction remains the last uncalculated two-electron 
second-order QED contribution for this state. 

This paper is organized as follows. In the next section we present the basic formulas 
for the two-photon exchange correction for the (ls)^2p3/2 state. The description of our 
numerical procedure is given in Sec. HI, and the results obtained are discussed in Sec. IV. 
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Relativistic units {h = c = 1) are used throughout this paper. 



II. BASIC FORMULAS 

The detailed derivation of the two-photon exchange corrections to the (ls)^2s and 
(ls)^2pi/2 states of Li-hke ions can be found in our previous paper p. For the (ls)^2p3/2 
state the derivation is performed along the same lines. Thus, we present mainly the final 
brmulas here. Our derivation is based on the two-time Green function (TTGFj method 
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16l |. For the detailed description of the method we refer to the recent review 17 1. 
The two-photon exchange corrections to the (ls)^2p3/2 state of the Li-like ions can be con- 
veniently separated in three parts: the two-photon exchange contribution due to the interac- 
tion between two Is electrons, the two-photon exchange contribution due to the interaction 
between the valence electron and one of the Is electrons, and the three-electron contribu- 
tion. The first part coincides with the two-photon exchange correction to the ground-state 
energy of He-like ions. Its calculation was carried out in [lO,lu|. This correction does not 
contribute to the 2p-2s splitting in Li-like ions and is not considered here. The remaining 
two-electron and three-electron corrections are diagrammatically depicted in Fig. ^ 

We start from the expression for the second-order correction to the energy shift of the 
level A; |l7|. 



A£^f = — 

" 27rz 



where Agkk{E) = gkk{E) - gfl{E), gkk{E) = {uk\g{E)\uk) , Uk is the unperturbed wave 
function, /S.E = E — E^\ E^^ is the unperturbed energy of the state k, and g^^^{E) = 
{AE)~^ is the function gkk{E) in the zeroth-order approximation. The function g{E) = 
g{E, x.[, ■ ■ ■ , x^; xi, ■ ■ ■ , x/yr) is the temporal Fourier transform of the N-electron two-time 
Green function. Its definition and the corresponding Feynman rules can be found in jl7|. 
The superscripts in Eq. (Q) indicate the order of the contribution in a. 

For the two-photon exchange correction, Feynman diagrams contributing to Ag^'^\E) 
are presented in Fig. ^ We refer to the corresponding contributions as the ladder (a), 
the crossed (b), and the three- electron (c) terms. The second term in Eq. is known as 
the disconnected contribution. It vanishes completely when considered together with the 



reducible contribution (for details, see [8|). In our case, the unperturbed wave function is 

Uk = ^(-l)^^/'Pa(xi)^p;,(x2)V'P,;(x3) , (2) 

Vol p 

where v denotes the valence electron, a and b are the electrons in (Is)^ shell, and P is the 
permutation operator (in the factor (—1)'^, the parity of the permutation is implied by P). 
For brevity we will use also the following notations: 

I{uj) = e^a'^a'D^^ito) , (3) 
labcdi^) = {ab\I{uj)\cd) , (4) 

Iab;cd = Iahcd{^bd) — hacd{^ad) , (5) 

r(.)^^, (a) 

where Aab = — ^fe, = (1, ck) are the Dirac matrices, and D^yioj) is the photon propa- 
gator. 

We separate the contributions of the diagrams under consideration into two parts: the 
reducible, with the energy of the intermediate state coinciding with the energy of the initial 
(final) state, and the irreducible, for the remainder, respectively. Omitting the derivation 
similar to that in Ref. Q], we present here only the final expressions for the energy shift. 
The reducible ("red") and irreducible ("ir") three-electron contributions read 

PQ 

E, Ip2P3nQ3{'Ap3Q3) I PlnQlQ2{AQlPl) /-\ 
n ^Ql + ^Q2 — ^Pl ~ 

^-^red = X! Kaav{^va){Iab;ab — hv;bv) 

(8) 



^^a 



~^ 2'^avvbi'^'va)Ibv;av + '^'^bvvai^va) -^va't^b 



where P and Q are the permutation operators, and the prime in the sum in Eq. ([7j) indicates 
that terms with the vanishing denominator should be omitted in the summation. In Eq. (jH} a 
and b denote Is electrons with opposite angular-momentum projections fia = —fib, v stands 
for the valence state with the angular-momentum projection fiy, and v is the valence state 
with fiy = 2/ia + /i^, (the corresponding contribution is assumed to be zero when is out of 
the range -j^, . . . , j„). 
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The irreducible two-electron contribution is 

AE^^, = AE^tr + Ai?^A + AE^ + AE-,^ , (9) 



riin2 



„,„„ Z/l J— oo 



(10) 



n\n2 



(e^ - U; - Sn^u) {Sc + UJ - Sn^u) 

7 /-co 

{ec-uj - en,u){e^ - u - Sn^u) 

A^:Lh = -E'^/ 

nin2 



(11) 



(12) 



(e^ - u; - £„iM)(£:^ - u; - Sn^u) 
Here we introduced the labels "lad" and "cr" for the ladder and the crossed diagram, and 
"dir" and "exch" for the direct and the exchange parts. The other notations are: 

F^^f{uj,nin2) = Icvmn2{^^)lnin2cv{uj) , (14) 

F^l{uj,nin2) = Y Icn2mv{'^)Imvcn2{'^) , (16) 
= Y Ivn2mv{^)Iniccn2{^ - ^vc) , (17) 

u = {1 — iO), and the prime on the sum indicates that some terms are excluded from the 
summation. First of all, we omit the reducible contribution, i.e. the terms for which the 
intermediate two-electron energy 6^ + £n2 equals the energy of the initial two-electron state 
Ey + Ec- Those are: (£„, Er,^) = (Sc^v) and (EyEc). In addition, we exclude also the infrared- 
divergent terms (see [a, for details), namely those with (e^je^J = (EcEv) in the direct 
crossed part and with (eni^n2) = i^c^c) and (EyE^) in the exchange crossed part. These terms 
should be considered together with the reducible contribution. Their sum can be shown to 
be infrared finite. We employ the notations Ai?>?fj,l and AE^^l^,, in order to emphasize that 
the corresponding terms are not "pure" irreducible and reducible contributions. 



We mention that the case under consideration differs from the cases of the 2s or 2ni/2 

n 

valence electrons considered previously in jSj by the fact that for the 2p3/2 Dirac state there 
is no adjoining state separated only by the finite-nuclear-size effect. Consequently, there is 
no need to exclude any further terms from the crossed contribution, as we had to proceed 
in Ref. in the case of the 2s and the 2pi/2 valence electron. 
Finally, we note the "reducible" contribution 



J-oo [UJ + lO)^ L 

+ cv) - + A.c cv) 

F\^^{uj - vc) - F\^^{-uj - A„„ vc) 

F!:^^{uj,vc)-F^^{-uj,vc)]. (18) 



III. NUMERICAL EVALUATION 

The three-electron contribution to the energy of (ls)^2s, (ls)^2pi and (ls)^2p3/2 levels 
of Li-like ions has been calculated in our recent investigation 3]. This evaluation appears as 
relatively simple since the corresponding expressions ((Zj), (jH)) contain at most one summation 
over the Dirac spectrum and no integrations over the virtual-photon energy. Thus we focus 
here on the calculation of the two-electron contribution. 

The summation over magnetic substates in Eqs. (fTUjl - ffTSj) . (fTHj) was performed by means 
of standard techniques. The resulting expressions can be found in Q]. As an independent 
check we employed also the direct numerical summation of Clebsch-Gordan coefficients. 

To calculate infinite summations over the spectrum of the Dirac equation in Eqs. (IIUII - 

n 

fll3|) . we employed the method of the B-spline basis set for the Dirac equation [20|. Typical 
basis sets contained 50 positive and 50 negative-energy eigenstates for each value of the 
angular-momentum quantum number k. The finite size of the nucleus has been taken into 
account employing the homogeneously-charged sphere model for the nuclear-charge distribu- 
tion. The values of the rms radii used in this work are the same as in B]. Infinite summations 
over n were truncated typically at \k\ = 10. Partial sums of the expansion over \k,\ were 
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fitted to the form 

= «o + E ^ (19) 

n=2 1^1 

using tlie least squares method. The coefficient yields the extrapolated value for the sum 
of the expansion. We found that different fits with N = 4-6 yield the same result with an 
accuracy of at least 5 digits. 

The integration over the energy of the virtual photon u in Eqs. fll0|) - ()13p represents the 
most difficult part of the calculation. To avoid strong oscillations for large values of u, we 
performed the Wick rotation of the integration contour. Deforming the contour, one should 
take care about the poles and the branch cuts of the integrand. The analytic structure of 
the integrand for Eqs. ()lH) - (fT!^ is shown in Figs. EE These graphs are very similar to 
those for the 2s- and 2pi/2-vdAence electrons in Ref. 0^. The only difference is that now 
three Dirac energy levels occur which are more deeply bound than the valence state: Is, 2s, 
and 2pi/2- The terms in Eqs. flll|) and (fTSj) containing these states and the valence state as 
intermediate were treated in a different way than the remainder, as is discussed below. 

For the evaluation of the direct parts of the ladder and crossed contributions, we perform 
the Wick rotation of the integration contours separating the corresponding pole contribu- 
tions, as shown in Figs. |21 and El In the direct part of the reducible contribution, we also 
perform a Wick rotation and then integrate by parts. This yields the following expression 
which can be evaluated directly, 

1 



^-^"red",dir — 2 



- - r du^l^^^2^P^ri^u;^ ' (20) 
n Jo + u'' duj 

where F'(A) = {dF/duj)^=A- 

Let us now turn to the exchange contribution. As one can see from Figs. |3]orEl in this 
case the integration contour is squeezed between two branch cuts of the photon propagators 
on the interval [0, A^d- Therefore, the standard Wick rotation of the contour is not possible. 
It is convenient to divide the contributions of Eqs. (jllj) and (jl3p into two parts. The first 
one accounts for the poles of the integrand on the interval [0, A^d and is referred to as the 
irregular part. The remainder is denoted as the regular part. This contribution does not 
possess any poles close to the squeezed part of the contour, which simplifies its numerical 
evaluation. However, it turns out as is the most time-consuming part of the calculation. 



One of the integration contours Creg used for the evaluation of the regular part is depicted 
in Fig. El 

The evaluation of the irregular part is less time consuming, but its structure is more 
difficult. In this case we need to take care of single and double poles of the integrand that 
are located close to the integration contour. The potential occurrences of one or two single 
poles and one double pole within the interval [0, A^d were treated by means of the following 
identities: 



UJ2 

duo — 

u)\ Xq 



uj±iO 



P 



-2 f{uj) . 

duj =1= t'KJ[Xo) 

u)i Xq — UJ 



duj 



Ll>2 



doj 



(xo — c<j ± iO)^ 



±z7r/'(a;o) + 



p 



Xq - UJi 
1 



(xq — ± iO)(xi — ± zO) 



Xi - Xq 
fU)2 



P 



(21) 



(22) 



Xq — U 



-P I duo 

J X\ — UJ 



=Fivr/(xo) ±mf{xi) 



(23) 



where P indicates the principal value of the integral. In Eq. the choice of the sign 

before i7r/(xo) and mf{xi) is determined by the sign of the infinitesimal addition ±iO in the 
first and the second denominator, respectively. For the numerical evaluation of the irregular 
contribution we employed the integration contour Cm shown in Fig. |31 It consists of 3 parts: 
[—ioo — e, — e], [— e, A„c + e], and [A^c + e, A^c + e + ^oo]. A small positive constant e was 
introduced in order to facilitate the numerical evaluation of the principal-value integrals. 
After integration by parts, the exchange contribution of the reducible part can be written 

as 



red" ,exch 



2m 



exch 



-P 



duj d \ 
-oo u duu 



per 



exch 



(At,C + UJ, cc) 



+FZ,^{u,vv)~2F'j^^{u;,vc) 



(24) 



It is again worth mentioning that the integral in Eq. (j24j) exists only if the sum of all 3 terms 
in the brackets is considered. For the each single term, the integral is infrared divergent. 
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IV. NUMERICAL RESULTS AND DISCUSSION 



The results of our calculations are presented in Table HJ where the direct, the exchange, 
and the three-electron contribution to the two-photon exchange correction of the valence 
2p3/2 electron with the (Is)^ shell are listed separately. The evaluation was performed 
within the Feynman gauge. We estimate the numerical uncertainty of our results to be 
less than 5 x 10~^ a.u. For bismuth, our results can be compared with the calculation by 
Sapirstein and Cheng j^. They report —6.529 and —6.670 eV for the two-electron and the 
three-electron contribution, respectively. This agrees well with our corresponding results of 
—6.5330 and —6.6698 eV, respectively. 

It is interesting to compare the results of the rigorous QED treatment with approx- 
imations evaluations based on relativistic many-body perturbation theory (MBPT). The 
difference between the QED and MBPT results can be conventionally regarded as a "non- 
trivial" QED contribution. In order to deduce the two-photon exchange correction within 
the framework of MBPT, we should introduce the following changes in our basic formulas: 
all summations over intermediate states should be restricted to positive-energy states only, 
the calculation should be performed within Coulomb gauge, and the virtual-photon energy in 
the photon propagator should be set equal to zero. Within this approximation, all reducible 
parts vanish, and the integration over the energy of the virtual photon can be carried out 
employing Cauchy's theorem. This yields zero for the crossed contribution, and finally we 
are left with the following expression for the total two-photon exchange correction within 
the MBPT approximation: 

AE^^Bprp = ' -^P2P3nQ3(0)-^PlnQlQ2(0) ^^5^ 

PQ e„>0 ^Ql + - £^P1 - £n 

(0) 

^-^MBPT — I _ ~ 5 V^OJ 

£„-^e„2>0 fcc + fci) fcni fcn2 

where the photon propagators should be taken in the Coulomb gauge and the prime indicates 
that terms with vanishing denominator should be omitted. We mention that Eqs. 
and (|26p include the contribution due to the exchange by two Breit photons (the B x B 
term). Strictly speaking, this term is of higher order than the level of validity of the Breit 
approximation, and, therefore, it appears to be inconsistent to include it within the MBPT 
scheme. 

In Table |n] and in Fig. IB] we compare the results of the rigorous QED treatment of the 



two-photon exchange correction to the 2p3/2-2s sphtting with the complete MBPT result 
[Eqs. (f^3j) and (j2SI)]; and the MBPT result dropping the B x B term. A similar analysis for 
the 2pi/2-2s splitting has been presented in our previous investigation [8]. Our results show 
that in the case under consideration the nontrivial QED contribution is essentially larger 
than that for the 2pi/2-2s transition. E.g., for uranium it yields —0.011 a.u., while the 
corresponding contribution to the 2pi/2-2s splitting is two times smaller, of about —0.006 
a.u. Moreover, we see that in our case the total correction changes its sign in the region 
between Z = 92 and 100. As a result, the MBPT result becomes incorrect by more than 
50% at very high values of Z. A further conclusion that can be drawn from our comparison 
is that for the 2p3/2-2s splitting the B x B term is of the same sign and magnitude as 
the nontrivial QED contribution. Thus, its inclusion improves the agreement between the 
MBPT and the QED result. This situation is contrary to the one for the 2pi/2-2s splitting, 
where the B x B term turns out to be of the same order of magnitude, but of different sign 
than the nontrivial QED contribution. 

To summarize this investigation we presented a rigorous QED evaluation of the two- 
photon exchange correction for the (Is) ^2^3/2 state of Li-like ions. Combining these results 
with the data for the (ls)^2s state from our previous study ^, we obtained the two-photon 
exchange correction for the 2p3/2-2s splitting. This is an important step towards the final 
goal consisting in the evaluation of all two-electron second-order QED corrections to the 
2p3/2-2s transition energy for the Li isoelectronic sequence. 
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TABLE I: Various contributions to the two-photon exchange correction for the (ls)^2p3/2 state 
of Li-Uke ions, in atomic units. The subscripts "dir" and "exch" label the direct and the exchange 
parts, respectively; the superscripts "2el" and "3el" refer to the two-electron and the three-electron 
contributions, respectively. 
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TABLE II: Comparison of the rigorous QED treatment of the two-photon exchange correction 
to the 2p3/2-2s splitting in Li-hke ions with the approximate MBPT treatment, in atomic units. 
B X B denotes the term corresponding to the exchange by two Breit photons. 
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FIG. 1: Feynman diagrams for the two-photon exchange corrections. 
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FIG. 2: The poles and the branch cuts of the integrand for the direct part of the ladder contribu- 
tion, and the integration contour Cld- 
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FIG. 3: The poles and the branch cuts of the integrand for the direct part of the crossed contri- 
bution, and the integration contour Ccd- 
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FIG. 4: The poles and the branch cuts of the integrand for the exchange part of the ladder 
contribution, and the integration contour Cjrr- 
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FIG. 5: The poles and the branch cuts of the integrand for the exchange part of the crossed 
contribution, and the integration contour Creg- 
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FIG. 6: The difference of the QED result for the two-photon exchange correction to the 2p3/2-2s 
transition and the corresponding MBPT results, with the B x B term included (solid line) and 
without this term (dashed line). The upper graph presents this difference in atomic units, and the 
lower one in units of per cent of the total QED contribution. 
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